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A  STUDY  OF  BLAST  LO.\DING 
TRANSMITTED  TO  BUILDING  FRAMES 

I .   INTRODUCTION 
A.   Intrcxiuctory  Statement 
The  loading  which  a  structure  experiences  due  to  the  pressxire  vaves  accom- 
panying blasts  frcm  high  explosive  and  atomic  bombs  has  been  the  subject  of 
much  study  and  research  since  the  use  of  such  bombs  was  begun  in  the  last  war. 
The  advent  of  more  powerful  atomic  bombs  and.   of  the  "super"  bombs  has  increased 
even  more  the  urgency  of  discovering  with  greater  exactness  the  nature  of  blast 
loads  on  structures. 

A  recent  study  in  the  field  of  blast  loading    has  concluded  that  a  gener- 
alized triangular  pulse  is  a  good  approximation  to  the  actual  transient  loads 
which  the  front  wall  of  a  structure  experiences.  However  depending  on  the 
action  of  the  wall  covering,  the  shape  of  the  load-time  curves  that  reach  the 
building  frame  itself  may  be  quite  different  frcm  the  pxilse  which  impinges  on 
the  wall  panels.   It  is  the  purpose  of  this  thesis  to  study  the  load-time  rela- 
tionship for  the  load  transmitted  by  the  wall  covering  to  the  supporting  frame. 
The  problem  is  necessarily  complex  and  is  approached  by  mailing  use  of  cer- 
tain simplifying  assumptions.  The  most  important  are: 

(l)  Given  a  typical  bay  of  a  one-story  frame-type  building,  consider  that 
the  wall  covering  of  the  bay  consists  of  a  "nest"  of  closely- spaced 
beams  of  unit  width  which  are  simply  supported  and  act  independently 
of  one  another.   (Figure  l)  These  beams  are  used  in  this  study  to 
approximate  the  action  of  a  slab  covering. 

^  ''n.  M.  Newmark,  "An  Engineering  Approach  to  Blast  Resistant  Design", 
Proceedings  ASCE,  Vol.  79,  Separate  No.  306,  October  1953. 
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(2)  Asetme  that  the  pressxire  locuiing  which  iiiQ)inge8  on  the  wall  panels 
can  be  represented  ty  a  triangular  pulse  function.   '  (Figure  2) 

(3)  AsBvaae  ideal  elas  to -plastic  action  of  the  beams.   (Figure  3) 

{k)  Assume  that  the  frame  members  supporting  the  wall  covering  remain 
fixed  in  space ^  at  least  during  the  inteirval  for  which  this  study 
is  made. 

It  is  doubtfvil  that  movement  of  the  frame  would  have  much  effect  on  the 
loads  reaching  the  frame  members.  A  deflection  of  the  frame  away  from  the 
blast  source  would  result  in  a  slight  decrease  in  the  beam  end  reactions; 
but,  since  the  x>eriod  of  vibration  of  the  supports  is  miich  larger  than  that 
of  the  wall  panels,  the  effect  should  be  negligible  except  possibly  for  load 
pulses  of  very  long  duration. 

If  the  beams  which  approximate  the  wall  covering  were  perfectly  rigid  and 
massless,  the  shape  of  the  loaxl  transmitted  to  the  frame  would  be  identical 
to  that  of  the  blast  pressure  pulse  (Figure  2).  However  in  our  case,  the  beams 
have  mass  and  elIso  have  elasto -plastic  characteristics;  therefore,  the  load 
which  reaches  the  supporting  members  is  much  altered. 

In  order  to  determine  the  frame  loading  for  the  given  assuniptions,  an 
elastic  analysis  and  a  simple  plastic  ajaalysis  are  \ised.  The  elastic  analy- 
sis is  derived  from  the  familiar  dynamic  beam  equation  (l)  and  the  plastic 
one  cooes  from  a  simplified  consideration  of  beam  action  when  a  plastic  hinge 
forms  at  the  center.  Prom  these  two  analyses  expressions  for  the  dynamic  end 
reaction  of  the  beams  on  the  siipports  are  derived. 

Frame  load  charts  giving  the  variation  of  dynamic  end  reaction  (frame  load 
per  unit  height)  with  respect  to  time  are  drawn  for  different  beam  resistance 


^  'N.  M.  Newmark,  "An  Engineering  Approach  to  Blast  Resistant  Design", 
Proceedings  ASCE,  Vol.  79,  Separate  No.  306,  October  1953- 
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characteristics  and  different  load  jnilse  diirations.  The  fraoe  loads  are  ccm- 
peired  graphically  with  the  load  which  the  beam  receives  as  a  resxilt  of  the 
hlast  pressure  wave.  Load -time  curves  for  the  supporting  frame  are  shown 
only  for  the  dvcration  of  the  loading  pulse;  however,  a  method  is  indicated 
t^  which  the  relation  may  be  found  for  times  beyond  the  end  of  the  prilse. 

As  brought  out  in  a  subsequent  section  of  this  study,  the  elastic  reaction 
expression  involves  the  sxmx  of  a  trigonometric  series.  For  certain  durations 
of  loading  pulses,  the  curves  for  one  term,  three  terms  euid  six  terms  are 
drawn  for  comparison.  By  means  of  this  graphical  comparison  of  a  few  repre- 
sentative pulse  durations,  it  is  demonstrated  that  the  use  of  only  the  first 
term  gives  results  which  axe  fairly  acc\irate  for  long  jmlses  but  not  very  eic- 
curate  for  short  xnilses.  For  all  the  remaining  pulse  durations,  only  the 
first  mode  (one  term)  is  xised  althoxjtgh  this  is  only  an  approximation  to  the 
exact  relation  which  theoretically  would  require  an  infinite  nimber  of  terms. 
No  attempt  is  made  by  analytical  means  to  evaluate  the  accuracy  obtainable  for 
various  pvilse  durations  by  using  only  the  first  mode.  Such  an  evaluation 
could  be  made  the  subject  of  further  study.  The  advantages  of  using  the  first 
mode  approximation  in  the  construction  of  the  frame  load  charts  is  apparent, 
however,  eind  such  use  is  intended. 

B.  Notation 

Terms  used  in  the  following  sections  are  defined  as  they  ajxpeax.  They 
are  assembled  here  for  convenience  in  reference.  Dimensions  of  the  quantities 
are  indicated  by  the  symbols  F,  M,  L,  T  in  parentheses  meaning  respectively 
force,  mass,  length  and  time. 

a^   =  El/m  (lW) 

9    =  angular  displacement 

9        "  angular  velocity  (1/T) 


K         =  1, 


0  =  angular  acceleration     (l/lT) 

r  ^ 

1  =  Impxilse  of  load  =  /    pdt   (FT/l) 
0(t)  =  time  distribution  factor 

length  of  beam  (wall  panel  strip)   (l) 

m    =  mass  per  unit  length  of  beam  (M) 

M_   =  dynamic  bending  moment  at  center  of  beam  (FL) 

M^   =  fully  plastic  moment  of  beam  (FL) 

M„   =  maximimi  static  center  moment  =  P  X    /&     (FL) 
S  m       ' 

Mj.   =  ultimate  moment  -  carrying  capacity  of  beam  (FL) 

p    =  instantaneous  value  of  applied  pressure  load  (f/l) 

p(x)  =  S  P  X  =  space  distribution  factor 
^^  n  n 

P    =  maximum  value  of  applied  pressure  load  (f/l) 
m 

R    =  dynamic  end  reaction  (F) 

R„   =  maximion  static  end  reaction  =  T  Jl /2     (F) 
5  m 

t    =  time  (T) 

t-        =  duration  of  load  pxilse     (t) 

t    =  time  of  beginning  of  second  elastic  phase   (T) 
e 

t    =  time  of  beginning  of  plastic  phase   (T) 

f   =  t-t   (T) 
p 

t"   =  t-t   (T) 
e 

T    =  natural  periods  of  vibration  of  beam  (t) 
n 

T^   =  fundamental  natural  period  of  vibration  of  beam  (T) 

T    =  a  function  of  time  alone 
n 

V    =  shear   (F) 

w    =  deflection  of  beam  (l) 


V  »  lineeir  velocity  of  beam  (L/T) 

w  =  distemce  from  center  of  mass  to  axis  of  rotation  (l) 

V  =  acceleration  of  mass  center  (l/T  ) 

w  «  effective  center  deflection  of  beam  at  onset  of  second 


e 


elastic  phase  (l) 


w    «  maximum  center  deflection  of  beam  which  occurs  at  t  «=  t   (l) 
c«  e  ^  ' 

m 

w    a  center  deflection  of  beam  at  t  »  t   (L) 
y  P 

X    »  characteristic  function  for  beam 

u    =  natural  (circular)  frequencies  of  beam  »  — s—  ^ —  (l/T) 

w,   =  fundamental  natural  frequency  of  beam  (i/T) 


II.   INITIAL  ELASTIC  BEHAVIC»  OF  WALL  PANELS 
A.  Statement  of  the  Problem 
1 .  The  Erynamic  Beam  Equation  ' 

VHien  a  beam  resists  moving  loads,  the  governing  differential  equation  of 
motion  for  the  beam  is,  according  to  the  ordinary  Bernoulli -Exiler  theory  of 
flexure , 

"^       0(t) 


Z^n  K 


2  Lw  _^  dfw  ^  a»  £iiLLiiiI= -2 (1) 


2 
where  a  «  El/ra;  w  =  deflection  measured  normal  to  the  beam  axis;  m  =  mass  per 

unit  length;  q  =  external  applied  unit  load;  x  =  distance  measured  from  origin 

at  left  end  of  beam;  and  t  =  time.  The  external  applied  load  is  defined  by  a 

m 
space  distribution  factor,  p(x)  «  \   P  X  and  a  time  distribution  factor,  p(t). 


^^^See  for  example,  S.  Tlmoshenko,  "Vibration  Problems  in  Engineering",  D.  Van 
Nostrand  Co.,  Inc.,  New  York,  N.  Y.,  Second  Edition,  July  1937,  p.  332. 


6 
Figure  4  is  a  simplified  sketch  of  the  loading  on  a  typical  beam.  For  the 
given  triangular  loading  pulse,  p(x)  =  1  since  the  beam  is  assumed  to  be  uni- 
formly loaded,  and  0(t)  =  P^(l  -  t/t J .   The  characteristic  function,  X  ,  for 

ml  n 

a  simply  supported  beam  is  sin  (nnx/^  ) .   The  constants,  P  ,  are  defined  as: 


i 


i 


p(x)  X  dx 


n 


Jo  "   ^ 


n 


Substituting  in  (2)  and  performing  the  integrations  leads  to: 

Pj,  =  ^  for  n  =  1,  3,  5,  7  


P  =  0  for  n  =  0;  2,  4,  6 


n 


The  solution  to  (l)  is  taken  in  the  form 


m 

w  =  )   X  X 

/    n  n 


(3) 


where  X  =  sin  (nnx/Ji   ),  a  function  of  x  alone,  and  where  x     is  a  function  of 
n  '         '  '  n 

time  alone.  The  substitution  of  (3)  in  (l)  with  the  subsequent  separation  of 
variables  leads  to  the  following  differential  equations : 

2  „ 

a 


2  ^n^^'^  ,,, 

T   +  U)    X   =  (5) 

n    n   n     m 


The  solutions  of  (h)   and  (5)  are  respectively: 


X  =  sin  —r-  (always  true  for  simply- 

•^^  supported  beam) 


7 
(6) 


t  «  A  cos  u  t  +  B  sin  u  t 
n    n      n     n      n 


Bl(j 


n 


sin 


0(t)  cos  u)  tdt  -  cos  (J  t  I   0(t)  sin  u  tdt 
'^        n  n  I   '^^  '      n 


(7) 


Equation  (7)  is  the  solution  to  the  single-degree  of  freedom  Equation  (^). 
The  first  two  terms  in  (7)  are  the  free  vibration  terms  and  the  last  two  the 
forced  vibration  terms.  For  initial  conditions  of  zero  displacement  and  zero 
velocity,  the  free  vibration  terms  disappear  since 


r 


w  (x)  X  dx 

o^  '  n 


*n  =  -T7 


¥0 


dx 


./. 


B  =  =- 
n   u 


J'   w  (x)  X  dx 
o   °     ° 


•'o 


dx 


where  w  (x)  euid  w  (x)  are  the  displacement  and  velocity  respectively  at  zero 
time. 


'^S.  Timoshenko,  "Vibration  Problems  in  Engineering",  D.  Van  Nostrand  Co.,  Inc., 
New  York,  N.  Y.,  Second  Edition,  July   1937,  p.  3^. 
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Substituting  Equations  (6)  and  (7)  into  the  expression  for  deflection  (3), 
yields  the  following: 


n   .  nnx 
w  =  )     sin  —r- 


n 


sin  u)  t  /  0(t)  cos  u  tdt 
n  /  '^^       n 

•'o 


cos 


w  t/" 

""J 
Jo 


0(t)  sin  ')  tdt 


(8) 


Substituting  the  values  of  P  and  0(t)  in  (8)  and  integrating,  results  in  the 
following  expression  for  w: 


4   1   m  .  nTOc 

V  = sm  —r- 

nn  mtj  (a)  t 
n  n 


n 


/  sin  u)  t 

(1  -  cos  >J  t)  -  r—  (1 .  .^ 

^                   n  '   t,  \  u)  t 

1  \  n 


Noting  that  w  for  a  simply- supported  beam  =  (n  n  / Ji    )  a.,   where  a  =^El/m: 


Ki      \      1  ,     nnx 

n=l,3,5--- 


(1  -  cos  u)^t)  -  :^  [  1  - 


sin  u)  t 
n_ 

(J  t 
n 


(9) 


B.  Reaction  and  Moment  Expressions 
1.   Pynamic  End  Reaction 

The  shear  at  any  point  x  along  the  beam  is  given  by: 


V  =  -  EI  ^ 


Performing  the  differentiation  on  (9)  and  evaluating  the  expression  at  the 
end  of  the  beam  (x  =  O),  the  dynamic  end  reaction  is: 


^  = 


m 


(l  -  cos  u)  t)  -  T—  [ 
°  ^1    \ 


n=l,3,5 — 


1  - 


sin  u  t 
n 

u  t 
n 


If  P   ,  the  fflfiLXlmum  value  of  the  external  load,  were  applied  as  a  static  load, 
the  maxinnjm  static  end  reaction  wo\ild  be: 

P  1 

For  ease  in  plotting,  the  reaction  expression  is  put  into  dimensionless  form 

2  2 

hy  noting  that  u  =  2n/T  =  2nn  /T,  since  T  =  T,/n  . 
n     '  n      '1       n    1 


sin  2nn  t/T,N 
l/(t  /T  )  .  t/T,  (  1 5 ^ 


(10) 


2,  Magnitude  of  Reaction  for  Long  Pulses 

As  the  duration  of  the  load  pulse  goes  to  infinity,  the  coefficient  of  the 
second  term  in  (10 )  approaches  zero  and  the  reaction  approaches: 


lim  =- 


-^  2  t 
cos  2nn  =r- 


n*l,3,5 — 


t^/T, 


The  maximum  reaction  then  approaches: 


10 


_„- •  '"7  /_n- ' '"''' ^"  i_7i- ^"7  ^'^^ 

ns=odd  n^odd  n«odd 

The  limit  of  the  summation  exists  as  shown.      It  is  given  hy  the  Rieznann  Zeta 
Function,   the  values  of  which  are  tabulated.      ' 
For  X  o  2,     C  (2)  =  TT^/6  =  1.645;   or: 


,.     ^       ,,     16     \      1        16 
t  /t,   — >•   •  n^odd 


(1  -  lA)  V 


2.00 


Then  for  pulses  of  duration  approaching  infinity,  the  dynamic  reaction  ap- 
proaches a  value  twice  the  maximian  static  reaction.  This  result  is  to  be 
expected  from  energy  considerations. 
3.  Magnitude  of  Reaction  for  Short  Pulses  and  Constant  Impulse 

For  short  pulses  of  duration  approaching  zero,  the  reaction  should  also 
approach  a  limit.  This  limit  can  be  found  from  a  consideration  of  the  reac- 
tion expression  which  holds  for  time  beyond  the  end  of  the  pulse  as  the  pulse 
diiration  becomes  very  small  and  total  impulse  remains  constant.  Figxore  5 
shows  the  new  time  origin  at  t  =  t. . 

As  before,  the  deflection,  w  =  ^   X  %   .     However  for  t  >  t. ,  the  forced 
'  '     ^/ n  n  1 

n 
vibration  terms  in  i^   drop  out,  leaving: 

T  =  A  cos  u  t  +  B  sin  u  t 

n    n  <     n     n      n 


^"''S.  Jahnke  and   F.  Emde,  "Tables  of  Functions  with  Formxilae  and  Curves", 
Dover  Publications,  New  York,  N.  Y.,  Fourth  Edition,  19^3. 
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The  initial  conditions  of  displacement  and  velocity  for  t  =  0  can  be  foxxnd 
from  Equation  (9)-   Setting  t  =  t  in  (9): 


/  \     m    \   1   .   nnx 

n=odd 


sin  u  t, 
n  1 

u  t^ 
n  1 


cos  u  t^ 
n  1 


Differentiating  (9)  with  respect  to  time  and  setting  t  =  t  ; 


vjx) 


Up  Jt 

m 


r? 


EI     /       n 


5  I 


cos  u   t,         , 
.  .  nil 

u     sin  u  t^    +  7 -  — 

n  n  1  t,  t. 


n=odd 


Values  of  A     and  B     are  then: 
n  n 


r^ 


\         V   (x)  X 
Jo         o'    '     n 


dx 


A     = 


"   n 


I: 


X  ^  dx 
n 


where  X     =  sin  —— 


A     = 
n 


UP  jL' 

n^EI     /       n^ 
n=odd 


sin  u  t,    -   u  t,   cos  u  t, 
n  1         n  1 n  1 

u  t^ 
n  1 


./ 


B     =  ~ 


r  w  (x)  X 

I  on 

1      Jo 


dx 


n  ~  cj       /•£ 

'     /  X 

^o        ^n 


dx 


12 


B 


n 


itesodd 


tj  t,  sin  u  t,  +  cos  u  t,  -  1 
n  1      n  1        n  1 

n  1 


Then  for  t  »  t  -  t, : 


v(x) 


i^P_i 


m 


IT^EI 


1  ^.  nnx 


cos  U)  t 
n 


J.  ,  sin  u  t-  -  (J  t,  cos  u  t- 
n  1 


n=< 


+  sin  J  t  I  — 


t-  sin  uj  t-  +  cos  u  t-  -  1 
1      n  1       n  1 

U)  t, 

n  1 


As  before,  find  the  shear  and  evaluate  at  x  «  0  in  order  to  get  the  dynamic  end 
reaction: 

CO 

^ —     )     -2    ^°^  ^n*  I ;7T J 

n     /   n  \         n  1  / 


+  sin  u  t 
n 


J,  ,  u  t-  sin  u)  t-  +  cos  u)  t-  -  1 
♦  (  n  1 n  1 n  1 

u  t, 
n  1 


(n) 


For  ary  pulse  duration,  the  ii^ulse  of  the  load  is: 


/ 


tl 


I  =  /   pdt  =  1/2  T^t^ 
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Now  let  t,  — >•     0  keeping  the  total  Impulse  constant.   P  then  Increases 
1  m 

rapidly,  and  in  terms  of  I  is  given  by: 

m   t 

Replacing  P  in  (ll)  by  its  value  in  terms  of  I  and  taking  t,  inside  the  sum- 
m  1 

mat ion: 


^  = 


6li 

2   , 
n   /   n 

n=odd 


cos  U)  t 
n 


sin  u  t^  -  tj  t,  cos  u  t, 
g  1    n  1 n_  1 

n  1 


+  sin  u)  t 

n 


„  ,  w  t^  sin  w  t,  +  cos  w  t^  -  1 

*  /  n  1 II  1 n  1 

n  1 


Nov  let  t  approach  zero: 


limR^ 
0 


bi 


1 


In  * 

n 


n=odd 


--^      \     -^  sin  u^t 
n         ^n 

n=odd 


Noting  again  that  u  =  2n/T  and  replacing  I  by  (l/2)  P  t  ; 


0 


sin  u  t 

n__ 


n=odd 


Take  the  first  mode  by  letting  n  =  1: 


11* 


S      1       ^1 


(12) 


The  foregoing  derivation  has  shovni  that  for  the  conditions  of  very  short 
pulse  duration  (say  when  t./T.  <  l/3)  and  constant  Impulse,  the  oaximuia  elastic 
reaction  approaches  a  value  equal  to  (8/tt)  t,/T-.  A  plot  of  reaction  for 
t  >  t,  in  which  the  load  pulse  satisfies  the  two  conditions  should  verify  (12). 
h.     Dynamic  Center  Moment  and  Beam  Resistance  Paraneter 

The  ntcment  at  any  point  x  along  the  beam  is  given  hy: 


M  »  -  EI 


n2 

o  w 


Performing  the  differentiation  on  (9),  evalixating  at  x  =    j^/2  auid  noting  that 

P   /   /8  is  the  mELXimum  static  center  moment  if  P    were  applied  as  a  static 
m*^  m 

loeui: 


ns 


-  v'^V'i)  •  */'i  (^  - 


sin  2nn     t/T^ 



2TTn     t/T^ 


(13) 


A  plot  of  kJvl   versus  the  time  parameter  t/T^  furnishes  a  convenient  means 
of  defining  the  limits  of  elastic  and  plastic  beam  action  for  different  pulse 
durations.  As  in  the  reaction  plots,  the  first  mode  (n  =  l)  in  (13)  is  used 
as  a  good  approximation.  Given  a  curve  of  Mjj/Mg  versus  t/T^,  the  fully  plastic 
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moment  of  the  beam  cross- section  may  be  considered  to  have  been  reaxihed  for 
different  veaues  of  M^.   When  a  given  KJvL^   is  said  to  equal  Mp/NL,  the  time 
corresponding  to  the  particular  value  of  M^M  is  the  time  of  hinge  formation, 
t  .  The  beam  resistance  parameter,.  Mp/M_  may  be  thotjght  of  as  the  ratio  of 
the  yield  value  of  the  beam's  moment  resistance  to  the  maximum  moment  caused 
by  the  static  application  of  P  . 

Ill   PIASTIC  BEHAVIOR  OF  WALL  PAMELS 
A.   Statement  of  the  Problem 
1.  Assumptions  in  the  Simplified  Plastic  Analysis 

In  order  to  derive  a  fairly  simple  expression  for  beam  reaction  when  the 
beam  ceases  to  behave  elastically,  the  following  assumptions  are  made: 

a.  A  plastic  hinge  forms  at  the  center  of  the  beam  when  the  dynamic  center 
moment,  M_^  reaches  a  value  equal  to  the  fully  plastic  moment  of  the  sec- 
tion, Mp.  This  time  is  designated  t  . 

b.  Neglect  the  elastic  bending  of  the  beam. 

c.  Assxime  that  deflections  are  small. 

d.  Assume  that  M^  remains  equal  to  Mp  until  the  beam  velocity  becomes 
eqiial  to  zero.  At  this  point,  M_^  becomes  less  than  Mp,  the  center 
hinge  disappesurs,  and  the  beam  behaves  elastically  again. 

Figure  6  represents  the  idealized  plastic  behavior  of  the  beam  conforming 
to  the  above  assumptions. 

Figure  7  is  a  free  body  diagram  of  half  of  the  beam  from  which  equations 
of  motion  may  be  written  as  follows: 


r 


F  .  =  M  w  where  w  =  ^  0;  w  =  -^  G 


w 


^T„.x,«  v^«  ,^.1 -l(|)^E 


/3 

"3r 


B.     Reaction  Egression 
Proa  Eqtiation  (15 )  solve  for  the  angular  acceleration: 


Substitute   (l6)  Into  Eqixatlon  {ik)  and  solve  for  R: 
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-g-  -R  =  -^  •  Tj-  © 


ilk) 


(15) 


Q 


2U 

.i3 


^-M, 


(16) 


p/   3Mp 
R  ■  *-g—  +  -^  vhere  p  =  function  of  time. 


(17) 


In  Equation  (17)>  P  is  the  external  load  acting  on  the  beam  at  tine  t'  measured 
after  the  formation  of  the  center  hinge.  Figure  8  represents  the  tiiae  origins 
for  the  different  phases  of  beam  action  which  occur  dviring  the  load  pulse. 
With 


P  =  P(t' )  »  P 


m 


^t*  +   t 


1  - 


and  t'  =  t  -  t  ,  and  noting  again  that  P  £/2   is  the  maximum  static  reaction: 
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%      1 
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(18) 


Differentiation  of  (l8)  with  respect  to  the  time  parameter,  t'/T,  gives  the 
resvilt  that  the  slope,  S,  of  the  plastic  reaction  curve  depends  only  on  the 
pulse  duration,  or: 

S  =  -  ^  l/it^/T:^)   =  -  ^  if  K  =  l/(t^T^)  (19) 


C.  Duration  of  Plastic  Behavior 


1.  Expression  for  t 


As  stated  in  the  assumptions  for  the  plastic  analysis,  the  center  hinge 
disappears  when  the  beam  velocity  beccanes  equal  to  zero.  The  time,  t  ,  at 
which  the  beam  stops  moving  down  is  therefore  the  end  of  the  plastic  phase 
and  the  beginning  of  the  second  elastic  phase. 

In  order  to  find  t  ,  an  expression  for  the  angular  velocity  of  the  beam 
can  be  derived  which  holds  during  the  plastic  phase  (t  <  t  <  t  ).  Equating 
this  angular  velocity  to  zero  and  solving  for  the  time  will  yield  the  equa- 
tion for  t  . 
e 

Differentiating  Eqxiation  (9)  with  respect  to  time  and  evalxiating  at 
X  =  ■J^/2   gives  the  linear  velocity  of  the  beam  center,  w  during  the  initial 
elastic  phase.  For  small  deflections,  the  average  angular  velocity  of  half 
of  the  beam  would  be  =  2/ J,    w  .  Performing  these  operations  and  setting 
t  =  t  will  give  the  equation  for  the  initial  angular  velocity  at  t'  =0 
(t  =  t  ).   (Figure  8)  Thus: 
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e*  = 


iiii 


8p  l^ 
m 

n  ^EI   /  n^ 

n=cxid 


n-1 


u  sin  u)  t  -  r—  (l  -  cos  u  t  ) 


n     n  p   t 


n  p' 


(20) 


The  required  angxilau:  velocity  may  now  be  found  "by  integrating  (l6)  which 
was  derived  from  the  equations  of  motion  for  the  pleustic  phase. 
Noting  that  p  in  (l6)  » 


p(f)  -  P^ 


f  +  t 


1  - 


3P 


©•  = 


m 


ai 


t*  +  t 


1  - 


Integrating  and  evaluating  the  constant  of  integration  hy  using  (20); 


'-^l--W-^ 


N 


n-1 


"5" 


(aJ  sin  u  t  -  r—  I  1  -  cos  u  t 
n     n  p   t,  \         n  p 


n«=odd 


J 


Slnplify  1^  letting  n  =  1  (first  mode)  and  note  that  u)      »  it  El/m/'    : 
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3P_ 


e'=-f(*' 


3TTUi^ 


Sin  w^t  -  --r-  (l 
1  p   <*),&, 


-  cos  U).  t  ) 
1  V 


(21) 


When  (21)  is  eqviated  to  zero,  the  resulting  eqxiation  is  a  quadratic  in  t' 
which  can  be  solved  "by  the  quadratic  formula  giving  the  folloving  resvilt  when 
put  in  dimensionless  form: 


t  t 

—    s  — «-    4- 

T  T 

■^1  ^1 


-.^f::^.!.:^-! 
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t     T 
1  ^1 
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3^  1 
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1-  cos 
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(22) 


IV.  SECOND  PHASE  OF  ELASTIC  WALL  PANEL  BEHAVIOR 
A.  Statement  of  the  Problem 
1,  Application  of  Dynamic  Beam  Eq\jation 

Equation  (22)  defines  the  time  at  which  the  beam  again  behaves  elastically. 
The  derivation  of  this  second  elastic  reaction  expression  follows  the  same 
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method  sis  used  on  the  Initial  one;  hovever,  one  free  vibration  term  is  added. 

When  t"  =  0  (t  =  t  ),  the  initial  velocity,  w  "  of  the  beam  is  zero,  but  the 

initial  displacement,  w  "  has  a  finite  value.  An  examination  of  the  formulae 

for  A  ajad  B  (Section  II)  shows  that  the  first  free  vibration  term  in  Equation 

(7)  has  a  value  other  than  zero. 

For  time  measured  from  t  (t"  =  t  -  t  ),  the  elastic  deflection  is  still 

e  e  ' 

of  the  form: 


m 
W"(x)  =   Y"X   T 

Z_  "  " 


where  X  =  sin  nnx/Z  and 
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1    B  A    cos  u)  t"  +  — 
n    n     n     mu> 


»t" 


sin 


y^t"  r   0(t")  cos  u^f'df 


.t" 


-  cos 


u)  t"  /   0(t")  sin  cj  f'dt" 
n   /    ^^   '     n 


(23) 


where  P  =  k/nv   when  n  =  odd.   (See  Equation  (2).) 

The  constant  A  is  defined  as 
n 


A  = 


"^;i 


Jo 


(x)  X  dx 


n 
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X  ^  dx 
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(2U) 
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^^  TT  y  ^^ 

Now  let  w^"(x)  =  w   sin  -r-  where  w   is  the  effective  center  deflection  at 
o       c      /       c 
e     -*■        e 

the  beginning  of  the  second  elastic  phase.  At  first  glance  it  appears  that  a 
good  approximation  of  w  "(x)  shotild  be  given  by 

w  "(x)  »  w   sin  ~ 
o  ^  '    c      / 

where  w   is  the  center  deflection  of  the  beam  in  the  plastic  phase  at  the 

in 
instant  the  beam  reaches  its  maximum  deflection  in  the  positive  direction  (in 

direction  of  applied  load)  and  begins  to  rebovmd  elastically.  However,  it  is 
api)arent  that  the  maximum  anqplitude  of  the  elastic  center  deflection  is  much 
less  than  the  corresponding  maximum  amplitude  of  the  plastic  center  deflection, 
Pigxire  9  illiistrates  roughly  this  difference.  The  solid  curve  gives  an  idea 
of  the  actual  variation  of  reaction  with  center  deflection  throughout  the  two 
elastic  phases  and  the  single  plastic  phase.  The  dotted  curve  is  eui  approxima- 
tion to  the  actual  behavior  which  suggests  a  simple  means  of  obtaining  a  close 

estimate  of  w   .As  shown  in  the  figure,  the  maxlmimi  plastic  center  deflec- 
c 
e  _ 

tion,  w   is  much  greater  than  w  :  and  w  ,  the  center  deflection  at  the  time 
'  c  c       y 

m  ®  _ 

of  hinge  formation  (t  )  is  approximately  equal  to  w 

^  ^e 

The  final  expression  and  best  approximation  to  be  used  for  w  "(x)  is  thus: 

M/       \  J  TTX 

w      (x)  =  w     sin  -rr 
o  ^    '         y  / 

The  value  of  w     may  be  foxmd  from  Equation  (9)  by  setting  x  =    £ /2  ajid 

t  =  t  .     The  result  is: 
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(25) 
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Nov  in  X  =  sin  —ji—  {2h) ,    let  n  =  1.   A  then  reduces  to; 
n       X  n 


A  = 


^o 


^  2  TTX  , 

=  T"  I   V  sin  -r-   dx 


=  w 


The  final  expression  for  A  is  thus  identical  to  w  and  is  given  by  (25). 

Go  back  to  the  expression  for  i^  ,  Equation  (23),  and  substitute  the  value 


n 


of  0(t") 


jZ)(f)  =  P 


m 


t"  +  t 


After  carrying  out  the  integrations  in  (23),  combining  it  with  X  and  simpli- 
fying, the  expression  for  elastic  deflection  for  t  <  t  <  t  ,  is : 


w"(x)   = 
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1       .      nrrx 
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1  - 


>      (26) 


B.  Reaction  Expression 
Having  derived  Equation  (26)  for  deflection,  the  reaction  expression  for 
the  second  elastic  phase  may  be  found  as  before  by  evaluating  the  shear  at 
X  =  0.   In  dimensionless  form  for  plotting,  and  using  the  first  mode  only,  the 


reaction  equation  for  t  <  t  <  t  ,  is: 
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^  "  FTt" 
r   1 


t"  +  t 


(27) 


V.   RESULTS  AND  DISCUSSION 
A.  Frame  Load  Charts 
1.  Construction  of  Charts 

In  the  previous  sections  of  this  paper,  expressions  for  the  beam  end  re- 
actions have  been  derived  for  the  elastic  and  plastic  phases  of  beam  behavior. 
The  reaction  is  expressed  as  a  ratio  of  the  dynamic  end  reaction  to  the  maxi- 
fflum  static  end  reaction,  (ILVr  );  the  time  parameter  is  written  as  the  ratio 
of  the  elapsed  time  to  the  fundamental  natural  period  of  the  beam  (t/T, );  and 
the  beam  resistance  parameter  is  expressed  as  the  ratio  of  the  fully  plastic 
moment  of  the  section  to  the  maximum  static  center  moment  (M^M  ).  For  con- 
venience, the  following  table  lists  the  various  phases  of  beam  action,  the 
time  range  of  each  phase  (Figure  8),  and  the  governing  reaction  equation: 
Phase  Time  Range  Equation 


Initial  Elastic 
Plastic 
Second  Elastic 


o  <  t  <  t 


t  <  t  <  t 
P       e 

t  <  t  <  t, 
e        1 


(10) 

(18) 
(27) 


Figvires  23,  27,  and  29  contain  plots  of  all  three  phases.  The  remaining 
reaction  curves  include  only  the  initial  elastic  and  the  plastic  phases.  On 


2k 
the  same  graph  with  each  reaction  plot  is  shown  the  external  load  which  im- 
pinges on  half  the  beam.  Thus  at  t  =  0,  the  load  on  half  the  beam  is  P  /^/2; 


m 


the  maximum  static  end  reaction  is  also  P  //2;  therefore  the  Initial  ratio 

m 

of  the  half -beam  load  to  the  maximxmi  static  reaction  is  equal  to  one.  The 
external  load  line  thus  steirts  at  unity  at  t  =  0,  and  drops  in  a  straight 
line  to  zero  at  t  =  t, .  This  is,  then,  a  graphical  representation  of  the  load 
which  arrives  on  half  of  the  beam,  and  what  this  half  beam  transmits  on  to 
the  frame  as  end  reaction. 

The  moment  cvirves  are  used  to  find  the  value  of  t  for  a  given  pulse 
duration  t  ,  and  resistance  parameter  (Mp/M  ) .   (Figure  l6)  The  duration 
of  the  initial  elastic  phase  is  thus  determined  graphically  from  the  plot  of 
M_/M„  versus  t/T^  by  allowing  VL   to  equal  Mp  as  Mp  clianges  over  a  range  of 
values.  This  brings  in  the  shape  and  physical  characteristics  of  the  beam 
since  the  fully  plastic  moment  varies  with  the  shape  and  material. 

The  duration  of  the  plastic  phase  Is  fixed  by  Eqxiation  (22).  The  loci  of 
these  times  as  fovurid  by  solving  (22)  form  a  smooth  curve  and  are  plotted  on  the 
reaction  graphs.   It  will  be  noted  that  for  a  given  piilse  duration,  the  loci  of 
Equation  (22)  indicate  the  expected  resxilts;  that  is,  that  beams  of  low  resist- 
ance (small  Mp/JL)  have  long  plastic  phases  sind  vice  versa. 

In  the  construction  of  the  reaction  curves,  it  was  found  that  a  discon- 
tinuity always  existed  between  the  elastic  and  plastic  phases.  This  "Jump" 
in  reaction  value  is  due  to  the  relatively  rough  assumptions  of  the  plastic 
analysis  used,  the  main  error  being  the  neglect  of  elastic  bending  of  the  beam. 
Since  the  source  of  the  discontinuity  was  known,  the  reaction  values  were  cor- 
rected by  the  amount  of  the  initial  reaction  error. 

Plots  of  the  second  elastic  phase  were  made  only  for  a  few  representative 
values  of  t^/T.  and  Mp/M  as  follows: 
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a.  As  an  example  of  a  short  pulse 
"b.  As  examples  of  long  p\ilses 


0.8 

1.0 

23 

5.0 

l.lf 

27 

5.0 

1.0 

27 

10.0 

l.U 

29 

10.0 

1.0 

29 

2.  Comparison  of  First  Mode  and  More  Exact  Reaction  Curves^  ' 

As  stated  in  the  introductory  statement,  no  attempt  is  made  in  this  study 
to  find  out  by  analytical  procedures  the  accuracy  with  which  the  first  mode 
approximation  represents  the  exact  elastic  reaction  curve.  However,  for  the 
initial  elastic  phase.  Equation  (lO)  was  expanded  and  data  was  collected  to 
plot  curves  for  one  term  (n  «=  l),  three  terms  (n  =  1,  3  and  ^),   and  6  terms 
(n  =  1,  3>  5^  7;  9  arid  11)  in  order  that  a  graphical  comparison  could  be  made. 

Figure  11  shows  the  results  of  the  plot  for  a  very  short  pvilse.   It  is 
obvious  that  the  first  mode  approximation  is  not  good  for  load  pulses  of  such 
short  duration.  Figure  35  illustrates  the  very  great  improvement  in  the  ac- 
curacy of  the  first  mode  approximation  for  ptilses  of  very  long  duration.  It 
woiold  seem  reasonable  to  expect  then,  that  as  the  load  pulse  duration  increases 
from  a  small  to  a  large  value,  the  accuracy  of  the  first  mode  should  show  a 
gradual  improvement.  Figure  19  furnishes  evidence  of  this  trend  since  it  shows 
the  first  mode  curve  to  be  a  better  approximation  than  it  is  for  t-i/T,  =0.1 
but  not  ajB  good  as  for  t,/T,  =  100.0. 


^    'This  subsection  refers  to  the  uncompensated  first  mode  reaction  cuirves  as 
derived  from  Equation  (10 )  by  setting  n  =  1.  These  curves  axe   indicated  in 
Figures  11,  19  and  35.  As  broiight  out  in  the  following  subsection,  aJLl  other 
reaction  plots  are   constrxjcted  -using  the  first  mode  multiplied  by  a  compen- 
sating factor. 
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One  situation  in  which  the  first  mode  curve  gives  a  false  impression  is 
encountered  at  very  small  values  to  t/T,  and  Kp/u^.     An  exaarple  is  Figure  11. 
Here  the  reaction  in  the  plastic  phase  for  Mp/M-  =  .05  appears  to  become  neg- 
ative before  the  end  of  the  pulse.  It  is  apparent,  however,  that  negative 
values  would  not  be  indicated  if  the  plastic  phase  plot  were  started  from  the 
more  exact  initial  elastic  curve  (for  3  or  6  terms).  An  examination  of  piilses 
of  longer  duration  will  indicate  this  same  ambiguity  at  very  small  values  of 
t/T^j^  and  Mp/Mg. 
3.  Method  Used  to  Make  First  Mode  Approximation  More  Conservative 

Figures  11,  19>  and  35  show  that  the  first  mode  reactions  plotted  frc«n 
Eqviation  (lO)  generally  furnish  val\ies  lower  than  those  indicated  by  the  more 
exact  cxurves.  In  order  to  retain  the  advantages  of  the  first  mode  approxima- 
tion but  at  the  same  time  to  obtain  more  conservative  end  reaction  magnitudes, 
the  reaction  expressions  for  the  two  elastic  phases  (Equations  (10 )  and  (27) 
with  n  =  l)  were  multiplied  by  a  factor  of  n  /8.  This  step  also  increased  the 
plastic  reaction  vaJ.ues  by  the  same  factor  since  the  plastic  phase  merges  with 
the  two  elastic  phases.  The  reasoning  behind  the  use  of  this  factor  follows. 

The  use  of  the  first  mode  curves  for  elastic  reaction  neglects  the  contri- 
bution of  the  higher  modes  in  (lO)  and  (27)  to  the  end  reaction.  Such  a  sim- 
plifying approximation  may  be  considered  as  being  equivalent  to  using  an  in- 
stantaneous beam  loading  represented  by  the  first  term  in  the  Fourier  Series 
expansion  of  a  uniform  load.  This  first  term  in  the  Fourier  expansion  is  merely 
a  simple  sine  curve  distribution  of  load.  Therefore,  the  actual  pressure  load 
which  at  any  instant  is  \jniformly  distributed,  is  replaced  by  a  sine  curve  load- 
ing the  total  ma-gnitude  of  which  is  very  close  to  that  of  the  actual  load. 
Figure  36  illxistrates  the  shapes  of  the  actual  ajid  approximately -equivalent 
loadings. 
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The  ratio  of  the  end  reaxition  calculated  from  the  correct  loading  to  that 
derived  from  the  equivalent  sine  curve  loading  is  found  as  follows:   (Figure  36) 

R 


(2)  =  2  J  n  P(^)  sin  -g-  dx  =  --til^ 

•'O  w 


then 


«(2)  ^ 


It  seems  reasonable  then,  that  reaction  values  more  nearly  equal  to  those 
resulting  from  a  uniform  load  will  result  if  Eqxiations  (lO)  and  (27)  (setting 
n  =  1)  are  multiplied  by  the  factor  n  /8.  Such  a  procedure  was  followed  in 
this  paper. 

Figures  11,  19  and  35  contain  plots  of  the  initial  elastic  reaction  expres- 
sion wherein  a  comparison  is  made  of  the  curves  plotted  from  Eqxiation  (lO)  for 
the  following  conditions: 

(1)  n  =  1 

(2)  n  =  1  using  multiplication  factor  of  n  /8 

(3)  three  terms  (n  =  1,  3,  5) 

(4)  six  terms  (n  =  1,  3,  5,  1,   9,  H) 

This  graphical  comparison  shows  that  the  use  of  the  first  mode  multiplied  by 
the  factor  n  /8  gives  results  which  are  more  conservative  approximations  to  the 
more  exact  curves. 

Since  veLlues  of  t  /T,  are  determined  graphically  from  the  ML/Mg  versus  t/T, 
plot,  another  soiu*ce  of  error  in  the  first  mode  reaction  curves  is  the  use  of 
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only  the  first  term  (n  =■  l)  in  Equation  (l3).  However,  the  elastic  center  mo- 
ment varies  as  — r-  for  n  equal  to  odd  integers;  therefore,  the  approximation  for 

n-^ 
center  moment  is  less  in  error  than  the  corresponding  approximation  for  reaction 

since  elastic  reaction  varies  as  -5-  .  A  conqjarison  of  the  center  moBients  due  to 

n 
the  two  loadings  in  Figiure  36  gives  the  following  result: 

"(1)  -^^^^-  0.125p(t)f2 


p 

M(2)  =  ^^^^l'^      (.500  -  .182)  =  0.129  p(t)/^ 
n 


then 


M/2)   .129   1.032 


.970 


By  the  foregoing  reasoning,  BK>re  accurate  values  of  t  /T,  could  be  gained 
"by  multiplying  Equation  (13)  "by  0.970.  This  factor  is  so  close  to  unity,  how- 
ever, that  little  change  in  the  values  of  t  /t,  would  result.  In  this  study. 
Equation  (13)  with  n  =  1  was  \ised  as  derived  without  applying  the  correction 
factor. 
k.     Discussion  of  Different  Phases  of  Wall  Panel  Action 

In  the  Introduction  to  this  paper,  the  statement  was  made  that  the  cheu^c- 
teristics  of  the  wall  x>anel8  would  cause  the  shape  of  the  load  treuasmitted  to 
the  frame  to  be  much  different  from  that  which  Impinges  on  the  pajiels  themselves, 
A  study  of  the  Frsuae  Loeui  Charts  (end  reaction  curves)  bears  out  this  statement. 
The  wall  panels  are  acted  vrpon  by  an  initial -peak,  triangular  loading  pulse 
but  send  to  the  frame  by  means  of  the  end  reactions,  a  loading  which  has  been 
radically  altered  in  shape. 
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The  unit  load  which  reaches  the  supporting  members  rises  gradually  from 
zero  along  the  initial  elastic  reaction  curve.   Depending  on  the  duration  of 
the  blast  pulse  and  the  beam  resistance  parameter,  the  beam  (wall  panel)  be- 
comes plastic  at  a  certain  time  and  the  reaction  then  decreases  linearly  with 
time  at  a  rate  depending  only  on  the  pulse  duration.  Then  if  the  peuael  has 
sufficient  elastic  resistance,  the  action  again  becomes  elastic  and  renains  so 
throvighout  the  remainder  of  the  load  pulse. 

For  very  short  pulse  durations,  the  plastic  phase  continues  to  the  end  of 
the  pulse  auad  no  second  elastic  phase  occurs.  Figures  11,  13^  and  15  illustrate 
this  behavior.  As  the  pulse  durations  increase,  the  second  elastic  phase  be- 
gins to  appear  starting  with  panels  of  high  elastic  resistance  (large  Mp/M^). 
Figures  27  and  29  show  how  (for  long  pulse  durations)  the  two  elastic  reaction 
curves  oscillate  about  the  incident  load  line.  It  should  be  noted  here  that 
damping  will  tend  to  decrease  the  an^jlitude  of  the  elastic  oscillations  about 
the  load  line.  If  the  effects  of  damping  could  have  been  included  in  the  de- 
rivations of  the  elastic  reaction  expressions,  Figures  27  and  29  would  proba- 
bly indicate  a  gradual  decrease  with  time  in  the  amplitude  of  the  vibrations  of 
elastic  reaction  about  the  external  load  line.  Thus,  in  the  second  elastic 
phase  for  low  beam  resistances,  the  elastic  reaction  would  become  coincident 
with  the  load  line  as  t  approaches  t  . 

Theoretically,  if  a  wall  panel  had  infinite  elastic  resisteuace,  its  end  re- 
action wovild  follow  the  initial  elastic  curve  and  the  frame  loeuis  would  rise 
to  a  maximum.  Plastic  action  takes  place  however,  since  the  panels  have 
finite  elastic  resistance,  and  prevents  the  frame  loads  from  building  up  to 
these  peeik.  values.  The  formation  of  the  center  plastic  hinge,  thus  effectively 
cuts  off  the  top  of  the  elastic  reaction  cxirve  and  consequently  permits  smaller 
loads  to  reach  the  frame. 
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The  occurrence  of  negative  values  of  reaction  in  the  elastic  oscillations 
about  the  load  line  may  be  highly  significant  under  certain  conditions.  If 
the  fasteners  used  to  attach  the  vail  panels  to  the  supporting  fraaie  are  not 
designed  to  withstand  tension  loads,  it  is  very  possible  that  such  fasteners 
will  fail  at  the  time  the  wall  panel  end  reaction  first  beccmes  negative. 
Such  failure  would  alter  the  frame  loaxls  in  the  case  of  load  pulses  of  suf- 
ficient duration  (t,/T-  >  0.8  in  the  frame  load  charts),  limiting  the  loads 
transmitted  by  the  jjanels  to  the  supporting  members  to  a  single  positive  pulse. 

B.  Tables  and  Graphs  of  Results 

Table  1  is  presented  for  convenient  reference.  For  given  values  of  pulse 
dxoration  (t,/T,  ),  plots  are  aJlso  made  of  the  beam  resistance  i)arameter  (Mp/M^,) 
versus  both  the  duration  of  the  initial  elastic  phase  (t  /T, )  and  the  duration 
of  the  plastic  phase  (t  /t,  -  t  /T, ).  (Figures  37  and  39)  Similar  graphs  are 
drawn  with  M_/m_  as  the  fixed  parameter.   (Figures  38  and  ko) 

Figure  37  shows  that  the  dui*ation  of  the  initial  elastic  phase  increases 
almost  linearly  with  increasing  values  of  Mp/M^  for  fixed  load  durations. 
Figure  38  demonstrates  that  for  a  given  beam  resistance,  the  duration  of  the 
load  pulse  has  little  effect  on  the  length  of  time  which  the  panel  remains 
elastic  except  in  the  case  of  the  shorter  load  pulses. 

As  might  be  exixjcted,  there  is  a  very  great  increase  in  the  duration  of  the 
plastic  phase  (for  given  pulse  duration)  with  a  decrease  in  waJJL  panel  resist- 
ance.  (Figure  39)  The  increase  is  most  pronounced  for  long  load  pulse  dura- 
tions. 

Figure  kO   shows  that  the  effect  of  load  duration  is  largest  for  low  values 
of  jjanel  resistance,  and  that,  in  general,  the  plastic  phase  duration  decreases 
as  the  dtaration  of  the  load  decreases  (with  given  panel  resistance). 
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C,  Frame  Loads   After  Termination  of  Load  Pulse 
For  time  beyond  the  end  of  the  pulse,  the  derivation  of  the  elastic  reaction 
expression  is  similar  in  procedure  to  that  used  in  Section  II  in  finding  the 
limiting  reeu:tion  for  short  pulses  and  constant  load  impulse.  Equation  (26) 
would  be  used  to  find  the  displacement  and  veloci-ty  at  t  »  t^  in  order  that 
A  and  B  could  be  evaliiated;  the  expression  for  w  would  then  be  found  where 


\      X  t 


and  X  still  equals  sin  rm-xj £  .  Reaction  could  be  found  as  before  by  evaluating 


n 
the  shear  at  x  =>  0. 


VI.  SUMMARy  AND  CONCLUSIONS 
A.  Applications  of  Frame  Load  Charts 

A  method  has  been  fo\md  for  determining  the  load -time  relation  for  blast 
loetd  transmitted  by  the  wall  panels  of  a  structure  to  the  supporting  members. 
The  frame  load  charts  may  be  used  for  three  possible  types  of  wall  panel  be- 
havior; a  purely  elastic  action,  an  elasto -plastic  action,  and  a  brittle  action. 

For  a  given  pulse  duration,  and  purely  elastic  behavior,  the  initial  elas- 
tic phase  reaction  curves  apply  throvighout  the  duration  of  the  pulse.  This 
assumes,  of  course,  that  the  wall  panel  has  sufficient  elastic  resistance  (high 
Hp/wu)  to  behave  in  this  manner. 

A  more  usual  wall  panel  behavior  is  the  elasto -plastic  I^t®'  Here  the  elas- 
tic resistance  is  not  high  enoiigh  to  confine  the  action  to  a  purely  elastic  one. 
Given  the  pulse  duration  ajid  beam  resistance  parameter,  the  frame  losul  charts 
xised  in  conjunction  with  the  dynamic  center  moment  curves  Indicate  the  loads 
which  are  transmitted  by  the  wall  panels  to  the  frame. 
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Wall  panels  constrvicted  of  brittle  material  vill  behave  elastically  up  to 
the  xiltimate  moaent  carrying  capewiity  of  the  panel  (KL)  and  then  fracture  sud- 
denly. The  time  of  fracture  could  be  determined  from  the  dynamic  center  moment 
curves  by  setting  Mj.  =  M^  over  a  range  of  values  of  M...  The  frame  load  charts 
wo\Ald  indicate  a  rise  in  load  along  the  initial  elastic  curve  up  to  the  point 
of  fracture  at  which  time  the  reaction  would  drop  to  zero. 

In  order  to  use  the  frame  loeui  charts  for  ^n  three  types  of  wall  panel  be- 
havior, the  following  information  must  be  available: 

a.  Weill  panel  span  and  a  curve  of  side-on  overpressure  versus  time 
from  which  P  ,   t. ,  R  ,  and  NL^  can  be  calcvilated .   ' 

b.  Dimensions  and  properties  of  wall  peuiel  in  order  that  T, ,  M- 
and  M.-  may  be  calculated  or  estimated. 

Entering  the  charts  with  t^/T,  and  Mp/KL  or  M.yM_,  the  magnitude  and  time 
distribution  of  the  uniform  frame  load  can  be  determined. 

B.  CONCLUSIONS 
The  frame  load  charts  presented  in  this  study  represent  an  attempt  to  eval- 
uate the  loads  received  by  the  supporting  members  from  the  building  covering 
when  the  walls  are  subjected  to  a  blast  pressure  wave.  Rough  assumptions  had 
to  be  made  in  some  cases  in  order  to  obtain  fairly  simple  expressions  from  which 
the  load  charts  could  be  constructed.  Also,  since  the  elastic  reawition  expres- 
sions involved  the  sum  of  an  infinite  series  of  terms,  they  could  not  be  used 
in  their  exact  form.  Therefore,  the  first  mode  approximations  to  the  exact 
elastic  reaction  expressions  were  modified  by  a  compensating  factor  of  tt  /8  and 
then  used  to  construct  the  fraiae  load  charts.  It  is  felt  that  this  step, 

TtT 

^''N.  M.  Newmark,  "An  Engineering  Approach  to  Blast  Resistant  Design", 

Proceedings  ASCE,  Vol.  79,  Separate  No.  306,  October  1953. 
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althoiigh  not  theoretically  exact,  provides  good  engineering  approximations  to 
the  actual  blast  loads  which  a  frame  would  experience.  Also,  keeping  in  mind 
the  additional  condition  that  uncertainties  in  loading  and  in  the  structural 
properties  of  the  wall  pejiels  will  always  exist,  the  author  feels  that  the 
frame  loads  indicated  hy  the  charts  are  sufficiently  accurate  to  "be  used  in  the 
design  of  building  frames. 

The  ba^ic  relations  for  panel  end  reaction  in  all  phases  of  elastic  and 
plastic  action  which  occur  diiring  the  load  pulse  have  been  derived.  For  the 
values  of  pulse  duration  and  beam  resistance  parameters  considered,  the  load 
charts  are  ccoplete  through  the  initial  elastic  and  plastic  phases  of  wall 
panel  action,  emd  representative  ciarves  for  the  second  elastic  phase  have  been 
drawn.  Thtis  the  main  purpose  of  the  study  --  that  of  determining  the  shape  of 
the  loads  transmitted  to  the  frame  —  has  been  accomplished.  Additional  work 
is  Inquired  to  finish  the  second  elastic  phase  and  to  extend  the  study  to  time 
beyond  the  end  of  the  pulse  using  the  methods  outlined  in  this  paper. 
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Figure  2-  Triangular  Load  Pulse 
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Figure  4  — Typical  Beam  with  Loading 
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Figure  5-  New  Time  Origin  for  t  >  t^ 
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Figure  6 -Plastic  Behavior  of  Beam 
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Figure  7  -  Free  Body  Diagram  of  Plastic  Beam 
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Figure  8  -  Time  Origins  for  Elastic  and  Plastic  Phases 
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Figure  9  -  Graphicea  Representation  of  Actual  and  Approxljnate  Variation 
of  End  Reaction  with  Center  Deflection 
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Figure  10  -  Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center 
Moment  Versus  Time  Parameter  (t^/T-,  ««  O.l) 
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Figure  11  -  Ratio  of  Dynamic  End  Reaction  to  Maximum  Static  End  Reaction  Versus  Time  Parameter 
(t^/T^  -  0.1) 
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Figure  12  -  Ratio  of  Dynamic  Center  Moment  to  Maximian  Static  Center 
Mcanent  Versus  Time  Parameter  (t./T,  =  0.2) 
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Figure  13  -  Ratio  of  Dynamic  End  Reaction  to  Maximum  Static  End  Reaction  Versus  Time  Parameter 
(tj^/Tj^  =  0.2) 
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Figure  lU  -  Ratio  of  Dynamic  Center  Mcment  to  Maximum 
Static  Center  Moment  Versus  Time  Parameter 
(t^/T^   -  0.3) 
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Figure  l6  -  Ratio  of  Dynamic  Center  Moment  to  Maximvmi  Static  Center 
Moment  Versus  Time  Parameter  (t^/T,  =  O.lf) 


^9 


50 


1.1* 


1.2  - 


1.0 


0.8  - 


Vm, 


0.6  - 


0.4  - 


0.2  - 


0.50 


Figure  18  -  Ratio  of  Dynamic  Center  Mament  to  Maximum  Static  Center  Moanent  Versus 
Time  Parameter  (t^^/T^  =  0.5) 
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Figure  20  -  Ratio  of  Dynanxic  Center  Moment  to  Maximum  Static  Center 
Moment  Veraiis  Time  Parameter  {\/'^2.   "  ^'^^ 
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Figure  21  -  Ratio  of  Pynamlc  End  Reaction  to  Maximum  Static  End  Reaction  Versus 
Time  Parameter  (t-,/T,  «  0.6) 
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Figure  22  -  Ratio  of  Dynamic  Center  Moment  to  Maalmum  Static  Center  Moment  Versus 
Time  Parameter  (t^/T^^  =  0.8) 
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Figure  2k  -   Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center  Mcoent  Versus 
Time  Parameter  (t^/T^  =  l-O) 


>l 


OJ 

O 

CO 

VD 

^ 

CVJ 

CVJ 

H 

H 

O 

O 

o 

o 

o 

V 

00 


> 


58 


2.0 


1.8  - 


1.6 


l.k 


1.2 


1.0 


V«S 


0.8 


0.6 


O.k 


0.2 


0.U8      0.56 


Figure  26  -  Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center  Moment  Versus  Tine  Parameter 
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Figure  28  -  Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center  Moment  VereuB 
Time  Parameter  l^^/'^i   =  10. 0) 
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Figure  30  -  Ratio  of  pynamic  Center  Moment  to  Maxlmvna  Static  Center  Moment  VerauB 
Time  Parameter  (t,/T,  »  20.0) 
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Figure  32  -  Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center  Moment  Versus 
Time  Parameter  (t,/T^  =  50. O) 
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Figure  3U  -  Ratio  of  Dynamic  Center  Moment  to  Maximum  Static  Center  Moment  Versus 
Time  Parameter  (t^^/Tj^  =  1CX).0) 
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Figure  36  -  Cooipaxison  of  Actual  and  Equivalent  Beam  Loading 
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Figure  37  -  Plot  of  IXrratlon  of  Initial  Elastic  Phase  (t  /T. )  Versus  Beaa  Resistance  Parameter  (M^M  ) 
for  Fixed  Values  of  Pulse  Duration  (tj^/T^^) 
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Figure  38  -  Plot  of  Duration  of  Initial  Elastic  Phase  (t  /t, J  Versus  Pulse  Duration  (t^/T, )  for 
Fixed  Values  of  Beam  Resistance  Parameter  (Mp/MI; 


71 


03 

VD 

^ 

CVJ 

o 

Ol 

W 

CVJ 

w 

OJ 

CO 


\£) 


CO 

d 


d 


o 

8 


S 


d 


o 


CVJ 

d 


I 

> 

I 


»4 


>. 


%  s 

"  I 

n 
«> 


£ 


V 


a 

•H 
4* 


^ 


72 


-n — I — r 


Mill  I  I — r 


?tn  I  I  r — I — r 


H 


I  I   I    I     I      I I 1 


'  ■  '    '    ■     ■      ' 


d 


-J    O 


o 

§ 


o 
d 


r- 


(0 

u 

V 

> 
a 


(0 

o 

a 
o 


o 
o 

H 


■■■'y'^mfi 


D  I  - 

FEB       2 

:am^  10  39 

1039 

MAR       I  6 

U  »♦  36 

MAY     20     • 
MY  2  1  56 

K'.Y3  1  6  1 

I  039 

10  5  9 

1                       Si^^V-          23306 

1  Thesis       f  ranc7 

1?V81             A  study  of  blast  load- 
1                   ing  transmitted  to  build- 
1                   ing  frames. 

1 

■              FEB       2 

DISPLA^I 
-^  —  10  3  9 

06^1 

■  UAK      16 

I                 MAY     20 

■  M  T  2  i  i  6 

n              MY  3  161 

10  39 

«*i*36 

1039 

10  3  9 

d 

1 

? 

rrancy 
yal                 A  study  of 

blast  loading 

25306 

tr?»r.§' 

fflitted  to  building  framas. 


